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ABSTRACT 

We  report  the  spontaneous  formation  of  a GaP  fullerene  cage  in  ab-initio  Molec- 
ular Dynamics  simulations  starting  from  a bulk  fragment.  A systematic  study  of  the 
geometric  and  electronic  properties  of  neutral  and  ionized  III-V  (GaP,  GaAs,  AlAs, 
A1P)  clusters  suggests  the  stability  of  hetcro-fullcrenes  formed  by  compounds  with 
zincblend  bulk  structure.  Our  prediction  is  supported  by  several  indicators:  these 
clusters  show  closed  electronic  shells  and  relatively  large  energy  gaps;  the  ratio  be- 
tween the  cohesive  energy  per  atom  in  the  cluster  and  in  the  bulk  is  very  close  to  the 
value  found  for  carbon  fullerenes  of  the  same  size;  the  clusters  are  thermally  stable 
up  to  a temperature  range  of  1500-2000  K and  they  do  not  dissociate  when  ionized. 

INTRODUCTION 

The  discovery  of  carbon  fullerenes  and  nanotubes  has  opened  a completely  new 
field  at  the  borderline  between  chemistry  and  physics  leading  to  many  new  phenom- 
ena and  applications.  Most  efforts  to  identify  fullerenes  based  on  other  elements 
have  focused  on  BN  which  is  the  most  similar  to  carbon  and  exists  in  nature  in  the 
hexagonal  (graphite-like)  structure[l-8].  However,  the  observed  (nested-)cages  and 
wires[l,2]  do  not  present  the  characteristic  pentagonal  rings  of  carbon  fullerenes. 

Here  we  explore  the  possibility  of  fullerene  cages  based  on  typical  semiconductors 
of  the  III-V  family,  like  GaAs,  GaP,  or  AlAs,  which  do  not  possess  a graphite-like  bulk 
structure.  We  have  recently  shown[9]  by  means  of  ab-initio  Car-Parrincllo  Molecular 
Dynarnics[10],  that  a small  GaP  bulk  fragment  spontaneously  organizes  in  a cage 
formed  by  a different,  number  of  atoms  of  the  two  elements  arranged  as  in  carbon 
fullerenes.  Here  we  extend  the  analysis  of  the  geometric  and  electronic  structure  of 
cages  with  the  stoichiometry  III3,-Va;±4  for  various  III-V  clusters  up  to  a total  of  52 
atoms. 

COMPUTATIONAL  DETAILS 

Our  results  are  obtained  by  means  of  the  Car-Parrinello  approach[10]  using  a 
Density  Functional  in  the  Generalized  Gradient  Approximation  proposed  by  Becke 
and  Pcrdow[ll,12].  This  approximation  reproduces  the  experimental  cohesive  energy 
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of  typical  bulk  semiconductors  within  5 valence  to  conduction  band  excitation  ener- 
gies^]. We  use  nonlocal  norm-conserving  first-principles  pseudopotentials[14]  and 
expand  the  single  particle  wavefunctions  on  a plane  wave  basis  set  with  a cut-off  of 
12  Rydberg.  We  use  a periodically  repeated  cubic  simulation  box  of  24  A side,  a 
size  found  large  enough  to  describe  isolated  clusters.  The  electronic  optimization  and 
structural  relaxation  have  been  performed  using  damped  second  order  dynamics  with 
electronic  mass  preconditioning  schemc[15].  The  symmetry  of  the  equilibrium  struc- 
ture is  not  biased  but  it  is  reached  spontaneously  during  the  geometry  optimization 
starting  from  the  corresponding  regular  polyhedron. 

RESULTS 

In  ref.  [9]  we  have  shown  that  a bulk  like  cluster  with  41  atoms  with  tetrahedral 
bulk-like  coordinations  spontaneously  evolves  to  a Gai6Pi2  bonded  cluster  with  12 
pentagons  and  4 hexagons  during  the  simulated  energy  minimization. 

The  observed  spontaneous  formation  of  a Ga]6Pi2  cage  with  pentagons  is  surpris- 
ing since,  in  the  case  of  BN[3-6],  (deformed)  squares  arc  found  to  be  energetically 
much  more  favorable.  For  B12Ni2[4],  there  is  an  energy  difference  of  9 cV  between 
the  cage  with  pentagons  and  the  one  with  squares  in  favor  of  the  latter  which  contains 
only  hcteropolar  bonds  and  is  favored  for  a material  composed  by  atoms  with  very 
different  electronegativity  as  B and  N.  Therefore,  most  studies  have  considered  cages 
B„,N„  formed  by  hexagons  closed  by  square  rings[3,4,6,7]. 

Very  recently,  Fowler  et  al. [8]  have  pointed  out  that,  among  the  cages  with  pen- 
tagons, those  with  one  species  in  excess  of  4 atoms  (B„N„+^)  minimize  the  number  of 
homopolar  bonds.  The  cage  GalcPi2  which  spontaneously  appear  in  our  simulation 
falls  into  this  class.  Therefore  we  focus  our  study  on  the  equilibrium  structure  and 
electronic  states  of  IIInVn±4  clusters,  with  a total  number  of  atoms  ranging  from  20 
to  52  and  for  different  III-V  compounds. 

In  Tables  I and  II  we  show  the  minimum  energy  structures  of  several  GaP  and 
GaAs  fullerenes  with  the  corresponding  symmetry  and  HOMO-LUMO  gap. 

As  in  the  case  of  carbon  and  BN  fullerenes,  the  III-V  clusters  would  represent 
metastable  states  with  respect  to  the  bulk  equilibrium  structure  and  only  experimen- 
tal observation  can  establish  with  certainty  their  existence.  We  support  our  prediction 
for  the  stability  of  the  examined  GaP  clusters  by  using  the  following  indicators:  i) 
closed  electronic  shells  and  large  energy  gaps;  ii)  cohesive  energy;  iii)  thermal  stability; 
iv)  stability  of  the  ionized  clusters. 

The  first  indicator  of  chemical  stability  is  the  energy  gap  between  the  highest 
occupied  and  lowest  unoccupied  molecular  orbitals  (HOMO  and  LUMO).  In  carbon 
fullerenes  a correlation  between  this  energy  and  the  observed  fullerenes  has  been 
experimentally  verified  [ 1 6] . 

In  figure  2 we  show  the  calculated  HOMO-LUMO  energy  gap  for  all  the  GaAs 
clusters  studied  as  a function  of  the  total  number  of  atoms  (the  nuclearity).  The 
clusters  with  a majority  of  atoms  of  type  V show  the  largest  HOMO-LUMO  gaps. 
In  general,  the  gap  grows  with  the  cluster  size,  with  a particularly  high  value  for  the 
cluster  with  52  atoms  as  found  also  for  the  analogous  BN  clustcr[8]. 
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Table  I.  Optimised  structures,  symmetry,  HOMO-LUMO  energy  gaps,  cohesive 
energies  and  Ga-P  binding  energies  for  several  GaP  dusters.  The  bulk  gap  is  given 
for  comparison.  The  annealing  temperature  is  also  given  for  the  first  two  structures. 


A comparison  of  the  binding  energies  per  atom  between  the  GaP  cages  and  the 
zincblend  bulk  phase  of  this  material  is  possible  only  for  the  clusters  with  the  same 
number  of  Ga  and  P atoms.  From  the  results  of  Table  I,  we  find  that  the  cohesive 
energies  per  atom  for  Gai2Pi2  with  squares,  Ga]2Pi2  with  pentagons,  and  Ga]0Pio 
are  about  10%  lower  than  in  the  bulk.  This  result  is  very  close  to  that  found  for  BN 
and  carbon  fullerenes  of  the  same  size[3]. 
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Table  II.  Optimised  geometries,  HOMO-LUMO  gaps  and  bond  lengths  for  some  of 
the  clusters  of  the  Ga2.AsI+4  seric. 


We  have  also  studied  the  thermal  stability  of  two  clusters  with  very  different 
energy  gaps,  namely  Ga^Ps  and  Ga8PI2  (see  Table  I).  For  both  clusters,  we  have 
performed  two  annealing  cycles  of  about  3 ps,  up  to  1500  K and  up  to  2000  K. 
The  system  is  heated  with  a rate  of  2xl015K/s,  then  equilibrated  for  one  ps  at  the 
highest  temperature,  and  finally  cooled  down  with  the  same  temperature  change 
rate.  For  Ga8P]2  no  bond  breaking  or  structural  rearrangements  occur  in  both  cycles 
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and  the  structure  comes  back  to  the  same  minimum  energy  configuration  when  the 
temperature  is  lowered.  This  is  also  the  case  for  Ga^Ps  in  the  annealing  up  to 
1500  K,  whereas  at  2000  I<  some  structural  rearrangement  takes  place  leading  to  a 
distorted  structure  with  higher  energy  when  cooled  down.  These  results  indicate  that 
the  thermal  stability  is  correlated  with  the  width  of  the  energy  gap. 


Figure  2.  HOMO-LUMO  gap  of 
the  GaAs  clusters  as  a function  of 
the  number  of  atoms  in  the  clus- 
ter. Dark  (red)  dots  correspond 
to  Ga^As^  clusters,  light  (green) 
dots  to  GaI+,iAsI  clusters.  The 
“magic”  clusters,  that  is,  the  more 
stable  of  the  serie,  have  been  cho- 
sen. 
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Mass  spectrometry  experiments  use  the  difference  in  mass-to-charge  ratio  of  ion- 
ized atoms  or  clusters  to  select  them.  Therefore  one  basic  requirement  for  the  possible 
detection  of  such  clusters  is  that  they  remain  stable  also  when  ionized.  We  have  in- 
vestigated the  stability  of  some  positively  ionized  clusters,  [GagPi2]+,  [Gai2P].6]+  and 
[Gai6Pi2]+-  The  electronic  structure  remains  almost  unaffected  and  degeneracies  are 
broken  by  negligible  amounts  in  the  order  of  hundredths  of  eV.  Only  minor  structural 
distortions  occur  upon  ionization.  In  particular,  the  six  equivalent  P-P  (Ga-Ga)  bond 
lengths  split  into  three  different  classes.  Remarkably,  during  a molecular  dynamics 
run  for  [Ga8Pi2]+  we  observe  a dynamical  exchange  between  these  three  classes  of 
bond  lengths  with  each  other.  This  effect  produces  features  in  the  low  frequency 
vibrational  spectrum  in  the  range  30—  120  cm"1  which  might  be  detected  by  infrared 
multiphoton  ionization  spectra[17]. 

CONCLUSIONS 

In  summary,  our  MD  results  strongly  suggest  that  small  fullerenes  formed  with  III- 
V semiconductor  compounds  could  be  stable,  since  they  have  high  symmetry,  closed 
electronic  shells,  and  large  HOMO-LUMO  energy  gaps.  These  clusters  arc  thermally 
stable  also  when  ionized. 

The  cages  with  the  stoichiometry  IIIX-V3;+,1,  which  minimizes  the  number  of  ho- 
mopolar  bonds,  are  particularly  stable  arid  show  the  largest  HOMO-LUMO  gaps. 

We  hope  that  this  work  will  stimulate  experimental  groups  to  widen  their  search 
for  hetero-fullerenes  also  to  III-V  semiconductor  compounds. 


W2.8.5 


ACKNOWLEDGMENTS 


FB  acknowledges  financial  support  from  the  PPM  (Prioriteitsprogramma  Materi- 
alenonder/oek)  project  number  96PPM001. 

REFERENCES 

1.  D.  Goldberg,  Y.  Bando,  0.  Stephan,  and  K.  Kurashima,  Appl.  Phys.  Lett.  73, 
2441  (1998). 

2.  P.  A.  Parilla,  A.  C.  Dillon,  K.  M.  .Jones,  G.  Riker,  D.L.  Schulz,  D.S.  Ginley  and 
M.J.  Heben,  Nature  397,  114  (1999). 

3.  G.  Seifert,  P.  W.  Fowler,  D.  Mitchell,  D.  Porezag,  and  TH.  Frauenheirn,  Chem 
Phys.  Lett.  268,  352  (1997). 

4.  F.  Jensen,  and  H.  Toftlund,  Chem  Pliys.  Lett.  201,  89  (1993). 

5.  X.  Blase,  A.  De  Vita,  J-.C.  Charlier,  and  R.  Car,  Phys.  Rev.  Lett.  80,  1666, 
(1998). 

6.  M-L.  Sun,  Z.  Slanina  and  S-.L.  Lee,  Chem  Phys.  Lett.  233,  279  (1995). 

7.  S.  S.  Alexandre,  M.  S.  C.  Mazzoni,  and  H.  Chacham,  Appl.  Phys.  Lett.  75,  61 
(1999). 

8.  P.  W.  Fowler,  K.  M.  Rogers,  G.  Seifert,  M.  Terrones,  and  H.  Terrones,  Chem 
Phys.  Lett  299,  359  (1999). 

9.  V.  Tozzini,  F.  Buda,  and  A.  Fasolino,  Phys.  Rev.  Lett.  85,  4554  (2000). 

10.  R.  Car  and  M.  Parrinello,  Phys.  Rev.  Lett.  55,  2471  (1985). 

11.  A.D.  Becko,  Phys.  Rev.  A 38,  3098  (1988). 

12.  J.P.  Perdew,  Phys.  Rev.  B 33,  8822  (1986). 

13.  G.  Ortiz,  Phys.  Rev.  B 45,  11328  (1992). 

14.  X.  Gonze,  R,.  Stumpf,  and  M.  Scheffier,  Phys.  Rev.  B 44,  8503  (1991). 

15.  F.  Tassone,  F.  Mauri,  and  R.  Car,  Pliys.  Rev.  B 50,  10561  (1994). 

16.  H.  Kietzmann,  R.  Rochow,  G.  Gant.efor,  W.  Eberhardt,  K.  Vietze,  G.  Seifert, 
and  P.W.  Fowler,  Phys.  Rev.  Lett.  81,  5378  (1998). 

17.  D.  Van  Heijnsbergen,  G.  von  Helden,  M.  A.  Duncan  A.  J.  A.  van  Roij  and  G. 
Meijer,  Phys.  Rev.  Lett.  83,  4983  (1999). 


W2.8.6 


